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We use Hubble parameter versus redshift data from Stern et al. [l| and Gaztanaga et al. 0] to place 
constraints on model parameters of constant and time-evolving dark energy cosmological models. 
These constraints are consistent with (through not as restrictive as) those derived from supernova 
Type la magnitude-redshift data. However, they are more restrictive than those derived from galaxy 
cluster angular diameter distance, and comparable with those from gamma-ray burst and lookback 
(3jT). time data. A joint analysis of the Hubble parameter data with more restrictive baryon acoustic 

oscillation peak length scale and supernova Type la apparent magnitude data favors a spatially-flat 
cosmological model currently dominated by a time-independent cosmological constant but does not 
exclude time- varying dark energy. 
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I. INTRODUCTION 
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It is well-established that the Universe is currently undergoing accelerated cosmological expansion. Observational 
evidence for the accelerated expansion comes from supernova Type la (SNIa) apparent magnitude measurements as 
a function of redshift [!, H[, cosmic microwave background (CMB) anisotropy data [|[ combined with low estimates 
of the cosmological mass density @, and baryon acoustic oscillation (BAO) peak length scale estimates [HE!!!!]- 
__i \ The underlying mechanism responsible for this accelerated expansion is not yet well characterized. The "standard" 
general relativistic model of cosmology has an energy budget that is currently dominated by far by dark energy, a 
negative-pressure substance that powers the accelerated expansion. (Another possibility is that the above observations 
are a manifestation of the breakdown of general relativity on large cosmological length scales. In this paper we assume 
that general relativity provides an adequate description of gravitation on cosmological length scales.) Dark energy 
^vjj , can vary weakly in space and evolve slowly in time, though current data are consistent with it being a cosmological 
■^j- ■ constant. For recent reviews see 

There are many dark energy models under discussion (for recent discussions see (Toj |. and references therein). 
The current "standard" model is the ACDM model [11| where the accelerated cosmological expansion is powered by 
Einstein's cosmological constant, A, a spatially homogeneous fluid with equation of state parameter wa = Pa/pa = — 1 
(where pa and pa are the fluid pressure and energy density). In this model the cosmological energy budget is 
dominated by far by p\, with cold dark matter (CDM) being the second largest contributor. The ACDM model 
provides a reasonable fit to most observational constraints, although the "standard" CDM structure formation model 
might be in some observational trouble (see, e.g., [l2j ). In addition, the ACDM model raises some puzzling conceptual 
H \ questions. 

If the dark energy density slowly decreased in time (rather than remaining constant like pa) , the energy densities 
of dark energy and nonrelativistic matter (CDM and baryons) would remain comparable for a longer period of time, 
and so alleviate what has become known as the ACDM coincidence puzzle. In addition, a slowly decreasing effective 
dark energy density, based on a more fundamental physics model that is applicable at an energy density scale much 
larger than an meV, could result in the current observed dark energy density scale of order an meV through gradual 
decrease over the long lifetime of the Universe, another unexplained feature in the context of the ACDM model. Thus 
a slowly decreasing dark energy density could resolve some of the puzzles of the ACDM model [l~3| . 

The XCDM parametrization is often used to describe a slowly decreasing dark energy density. In this parametriza- 
tion the dark energy is modeled as a spatially homogenous (X) fluid with an equation of state parameter wx = Px/px, 
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where wx < —1/3 is an arbitrary constant and px and px are the pressure and energy density of the A-fluid. When 
wx = — 1, the XCDM parametrization reduces to the complete and consistent ACDM model. For any other value of 
wx(< —1/3) the XCDM parametrization is incomplete as it cannot describe spatial inhomogeneities (see, e.g., fl3]). 
Here we study the XCDM parametrization only in the spatially-flat cosmological case. 

The 0CDM model — in which dark energy is modelled as a scalar field 4> with a gradually decreasing (in <$>) potential 
energy density V{4>) — is the simplest complete and consistent model of a slowly decreasing (in time) dark energy 
density. Here we focus on an inverse power-law potential energy density V(<f>) oc 4>~ a , where a is a nonnegative 
constant [13J, |l5|. When a = the 0CDM model reduces to the corresponding ACDM case. Here we only consider 
the spatially-flat 0CDM cosmological model. 

It has been known for some time that a spatially-flat ACDM model with current energy budget dominated by 
a constant A is largely consistent with most observational constraints (see, e.g., [1(1 H3)- SNela, CMB, and BAO 
measurements mentioned above indicate that we live in a spatially-flat ACDM model with nonrelativistic matter 
contributing a little less than 30 % of the current cosmological energy budget, with the remaining slightly more than 
70 % contributed by a cosmological constant. These three sets of data carry by far the most weight when determining- 
constraints on models and cosmological parameters. 

Future data from space missions will significantly tighten the constraints (see, e.g., |18j). However, at present, it is 
important to determine independent constraints that can be derived from other presently available data sets. While 
these data are not yet as constraining as the SNela, CMB and BAO data, they potentially can reassure us (if they 
provide constraints consistent with those from the better known data) , or if the two sets of constraints are inconsistent 
this might lead to the discovery of hidden systematic errors or rule out the cosmological model under consideration. 

Other data that have been used to constrain cosmological parameters include galaxy cluster gas mass fraction (e.g., 
[I3,|T9]), gamma-ray burst luminosity distance (e.g., [2(| El), large-scale structure (e.g., [HI): strong gravitational 
lensing (e.g., [13), and angular size (e.g., [1, HH [H[) data. While the constraints from these data are less restrictive 
than those derived from the SNela, CMB and BAO data, both types of data result in largely compatible constraints 
that generally support a currently accelerating cosmological expansion. This gives us confidence that the broad 
outlines of the "standard" cosmological model are now in place. 

Measurements of the Hubble parameter as a function of redshift, H(z) , have also been used to constrain cosmological 
parameters (see [26| for a review). A variant of this test uses lookback time data (see, e.g., [13, [HI). Building on 
the work of [2{J, Simon et al. [30] used the differential ages of 32 passively evolving galaxies to determine 9 H(z) 
measurements in the redshift range 0.09 < z < 1.75. Cosmological constraints derived using these data are described 
in [3l|, HH ; more recent references may be traced through (33| . 

Stern et al. (2010, hereafter S10) [IJ extended the Simon et al. 30] sample to 11 measurements of H(z) in the 
redshift range 0.1 < z < 1.75. These data have been used for cosmological tests by Shafieloo & Clarkson (34|. It 
has become common to augment the S10 data with the Gaztahaga et al. (2009, hereafter G09) Q estimates of H(z) 
determined from line-of-sight BAO peak position observations. These data, listed in Table 1, have also been used to 
constrain cosmological parameters (see, e.g., [H, [36[ and references therein). There are problems with a number of 
these analyses. Some of them include both the G09 data points at z = 0.24 and z = 0.43 (which we use here), as well 
as the G09 single summary data point at z = 0.34 that is based on exactly the same data as the two individual points. 
In addition, a number of these analyses either ignore the G09 systematic errors or incorrectly account for them. We 
account for the G09 statistical and systematic errors by combining them in quadrature; the G09 data points we list 
in Table 1 are identical to those used by Ma & Zhang [35( and Zhang et al. 26] . 

In this paper we use the 13 S10 and G09 H(z) measurements listed in Table 1 to constrain the ACDM and 0CDM 
models and the XCDM parametrization. The resulting constraints are compatible with those derived using other 
techniques. We also use these H{z) data in combination with BAO and SNela measurements to jointly constrain 
cosmological parameters in these models. Adding the H{z) data tightens the constraints, somewhat significantly in 
some parts of parameter space for some of the models we study. 

Our paper is organized as follows. In Sec. [TT] we present the basic equations of the three dark energy models we 
study. Constraints from the H{z) data are derived in Sec. IIIII In Sec. IIVI we determine joint constraints on the dark 
energy parameters from a combination of data sets. We summarize our main conclusions in Sec. |Vj 

II. BASIC EQUATIONS OF THE DARK ENERGY MODELS 

The Friedmann equation of the ACDM model with spatial curvature can be written as 

h 2 (z, h , p) = Hi [a m0 (i + Z ) 3 + n A + (i - n m0 - n A )(i + zf] , (i) 

where z is the redshift, H(z, Ho, p) is the Hubble parameter, Ho is the Hubble constant, and the model-parameter set 
is p = (fl m o,flA) where Sl m o is the nonrelativistic (baryonic and cold dark) matter density parameter and Q\ that 
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of the cosmological constant. Throughout, the subscript denotes the value of a quantity today. In this paper, the 
subscripts A, X and <fi represent the corresponding quantities of the dark energy component in the ACDM, XCDM 
and 0CDM scenarios. 

In this work, for computational simplicity, spatial curvature is set to zero in the XCDM and 0CDM cases. Then 
the Friedmann equation for the XCDM parametrization is 

H 2 (z, Ho, p) = H 2 [n m0 (l + zf + (1 - O m0 )(l + z) 3 ( 1+tu -)] , (2) 

where the model-parameter set is p = (fi m o, wx)- 

In the 0CDM model, the inverse power law potential energy density of the scalar held adopted in this paper is 
V(4>) — Knip(j)~ a , where m p is the Planck mass, and a and k are non- negative constants [l5j . In the spatially-flat 
case the Friedmann equation of the </>CDM model is 

Sir 

H 2 (z,H ,p) = —^(p m + P4> ). (3) 

Here H(z) = a/a is the Hubble parameter where a(t) is the cosmological scale factor and an overdot denotes a time 
derivative. The energy densities of the matter and the scalar field are 

m 1 i 

Pm = ^a~\ (4) 

and 

2 

771 

/^ = ^W> 2 + «m^-«), (5) 

respectively. According to the definition of the dimensionless density parameter, one has 

8irp m p m 

^ = 3^H 2 = p—^pT, (6) 

The scalar field 4* obeys the differential equation 

+ 3^-^m>-^ +1 >=O. (7) 
a 2 F 

Using Eqs. ^ and ([7]), as well as the initial conditions described in [l5|], one can numerically compute the Hubble 
parameter H(z). In this case the model-parameter set is p = (Q m o,a). 



III. CONSTRAINTS FROM THE H(Z) DATA 

We use the 13 H(z) measurements of S10 and G09 listed in Table 1 to constrain cosmological parameters. We 
constrain cosmological parameters by minimizing \h, 

Xhk h o,p) - 2 • W 

i=l °H,i 

Here Zi is the redshift at which H{zi) has been measured, H th is the predicted value of H(z) in the cosmological model 
under consideration and H ohs is the measured value. From Xn(Ho,p) we compute the likelihood function L(Ho,p). 
We then treat Hq as a nuisance parameter and marginalize over it using a gaussian prior with Hq = 68 ± 3.5 km s" 1 
Mpc -1 [13, HI] to get a likelihood function L(p) that is a function of only the cosmological parameters of interest. 
The best-fit parameter values p* are those that maximize this likelihood function and the 1, 2 and 3 a constraint 
contours are the set of cosmological parameters (centered on p*) that enclose 68.27, 95.45 and 99.73 %, respectively, 
of the probability under the likelihood function. 

Figures[T][3]show the constraints from the H (z) data on the three dark energy models we consider. Not unexpectedly, 
these show that the H(z) data fairly tightly restrict one combination of the cosmological parameters while leaving 
the "orthogonal" combination relatively unconstrained. Comparing these results to those shown in Figs. 1-3 of [32j, 
determined using the 9 H(z) data points of Simon et al. 30] , we see that the newer S10 and G09 data result in 
significantly more restrictive constraints on cosmological parameters. These constrains are comparable with those 
from gamma-ray burst (|ll|, Figs. 1-3) and lookback time (p7|. Figs. 1-3) data. They are more restrictive than those 
that follow from galaxy cluster angular diameter distance data ([2J], Figs. 1-3). 
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IV. JOINT CONSTRAINTS 

Following [24j . we derive constraints on cosmological parameters of the three models from a joint analysis of the 
H{z) data with the BAO data [3] and Union2 compilation of 557 SNela apparent magnitude measurements (covering 
a redshift range 0.015 < z < 1.4) from ||. 

Figures 2H6] show the constraints on the cosmological parameters for the ACDM and (fiGDM models and the XCDM 
parametrization from a joint analysis of the BAO and SNela data, as well as from a joint analysis of the BAO, SNela 
and H{z) data. Adding the H(z) data tightens up the constraints somewhat, most significantly in the ACDM case 
(Fig. 4) and least so for the </>CDM model (Fig. 6). 

Figures [7H9] display the one-dimensional marginalized distribution probabilities of the cosmological parameters for 
the three cosmological models considered in this work, derived from a joint analysis of the BAO and SNela data, as 
well as from a joint analysis of the BAO, SNela and H(z) data. The marginalized 2 a intervals of the cosmological 
parameters are presented in Table [TTJ 

The combination of BAO and SNela data gives tight constraints on the cosmological parameters. Adding the 
currently-available H (z) data to the mix does shift the constraint contours, however the effect is significant only in 
some parts of parameter space for only some of the models we study. While useful, current H(z) data do not have 
enough weight to significantly affect the combined BAO and SNela results in most parts of model-parameter space. 
The H(z) data have a little more weight than currently-available gamma-ray burst luminosity measurements ([HI, 
Figs. 4-6 and 10-12). 

In summary, the H(z) data considered here are very consistent with the predictions of a spatially-flat cosmological 
model with energy budget dominated by a time-independent cosmological constant. However, the data do not rule 
out time-evolving dark energy, although they do require that it not vary rapidly. 

V. CONCLUSION 

We have shown that the Hubble parameter versus redshift data from S10 and G09 can give interestingly restrictive 
constraints on cosmological parameters. The resulting constraints are compatible with those derived from other 
current data, thus strengthening support for the current "standard" cosmological model. The H (z) data constraints 
are approximately as restrictive as those that follow from currently-available gamma-ray burst luminosity data and 
lookback time observations, and more restrictive than those from currently- available galaxy cluster angular size data. 
They are, however, much less restrictive than those that follow from a combined analysis of BAO peak length scale 
and SNela apparent magnitude data. 

The spatially-flat ACDM model, currently dominated by a constant cosmological constant, provides a good fit to 
the data we have studied here. However, these data do not rule out a time-evolving dark energy. 

As discussed in [35j, future high-z, high-accuracy H(z) determinations from BAO observations will likely result in 
cosmological parameter constraints comparable to those that follow from SNela data. 
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FIG. 2: 1, 2, and 3 a constraint contours for the XCDM parametrization from the H(z) data. The dashed horizontal line 
at ujx = —1 corresponds to spatially-flat ACDM models. The star marks the best-fit pair (f2 m o,tox) = (0.29,-0.88) with 

Xmin = 10.1. 




FIG. 3: I, 2, and 3 a constraint contours for the <j)CDM model from the H(z) data. The horizontal axis at a = corresponds 
to spatially-flat ACDM models. The star marks the best-fit pair (n m o, a) = (0.28, 0.46) with Xmin = 10.1. 
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FIG. 4: Thick (thin) solid lines are 1, 2, and 3 a constraint contours for the ACDM model from a joint analysis of the BAO and 
SNela (with systematic errors) data, with (and without) the H(z) data. The cross ("+") marks the best-fit point determined 
from the joint sample without the H(z) data at fi m o = 0.28 and Qa = 0.76 with Xmin — 531. The star ("*") marks the best-fit 
point determined from the joint sample with the H(z) data at fl m o = 0.28 and Qa = 0.72 with Xmin = 541. The dashed sloping 
line corresponds to spatially-flat models. 
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0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 
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FIG. 5: Thick (thin) solid lines are 1, 2, and 3 a constraint contours for the XCDM parametrization from a joint analysis of 
the BAO and SNela (with systematic errors) data, with (and without) the H(z) data. The cross ("+") marks the best-fit point 
determined from the joint sample without the H(z) data at Q m o = 0.28 and lux = —1.04 with Xmin = 531. The star ("*") 
marks the best-fit point determined from the joint sample with the H(z) data at fl m o = 0.28 and u)x = —1.01 with Xmin = 541. 
The dashed horizontal line at lux = — 1 corresponds to spatially- flat ACDM models. 
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FIG. 6: Thick (thin) solid lines are 1, 2, and 3 a constraint contours for the 0CDM model from a joint analysis of the BAO and 
SNela (with systematic errors) data, with (and without) the H(z) data. The cross ("x") marks the best-fit point determined 
from the joint sample without the H(z) data at fl m o = 0.28 and a = with Xmin — 531. The diamond ("0") marks the best-fit 
point determined from the joint sample with the H(z) data at fl m o = 0.28 and a = 0.0 with Xmin = 541. The a — horizontal 
axis corresponds to spatially-flat ACDM models. 
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FIG. 7: One-dimensional marginalized distribution probabilities of the cosmological parameters for the LCDM model. Thick 
(thin) lines are results from a joint analysis of the BAO and SNela (with systematic errors) data, with (and without) the H(z) 
data. 
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FIG. 8: One-dimensional marginalized distribution probabilities of the cosmological parameters for the XCDM parametrization. 
Thick (thin) lines are results from a joint analysis of the BAO and SNela (with systematic errors) data, with (and without) 
the H(z) data. 
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FIG. 9: One-dimensional marginalized distribution probabilities of the cosmological parameters for the 0CDM model. Thick 
(thin) lines are results from a joint analysis of the BAO and SNela (with systematic errors) data, with (and without) the H(z) 
data. 
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TABLE I: Hubble parameter versus redshift data from S10 and G09. Where H(z) and an are in km s 1 Mpc 



Model BAO + SNela H{z) + BAO + SNela 

ACDM 0.24 < O m0 < 0.33 0.24 < n m0 < 0.33 

0.5 < Q A < 0.97 0.51 < Q A < 0.9 

XCDM 0.24 < O m0 < 0.33 0.24 < Q m0 < 0.33 

-1.30 < uj x < -0.80 -1.26 < lux < -0-79 

cj>CDM 0.24 < fi m0 < 0.33 0.24 < Q m0 < 0.33 

< a < 0.73 < a < 0.76 



TABLE II: Two standard deviation bounds on cosmological parameters. 



